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Abstract 

The influence of the temperature and its fluctuations on the ion saturation current and the 
floating potential, which are typical quantities measured by Langmuir probes in the turbulent edge 
region of fusion plasmas, is analysed by global nonlinear gyrofluid simulations for two exemplary 
parameter regimes. The numerical simulation facilitates a direct access to densities, temperatures 
and the plasma potential at different radial positions around the separatrix. This allows a compar- 
ison between raw data and the calculated ion saturation current and floating potential within the 
simulation. Calculations of the fluctuation-induced radial particle flux and its statistical properties 
reveal significant differences to the actual values at all radial positions of the simulation domain, 
if the floating potential and the temperature averaged density inferred from the ion saturation 
current is used. 



This is the preprint version of a manuscript submitted to Plasma Physics and Controlled Fusion. 
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I. INTRODUCTION 



In the edge and scrape-off layer regions of magnetically confined plasmas the fluctuating 
plasma density n = ric — n^, the plasma potential $ and the radial particle flux Fr are 
usually inferred from Langmuir probe measurements. However, the quantities measured 
by conventional cold Langmuir probes are the ion saturation current lis and the floating 
potential Vg. Following the elementary Langmuir probe theory, these are related to the 



density anc 
Te and Tj 



theplasma potential by expressions involving the electron and ion temperatures 
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A Maxwellian electron velocity distribution is assumed, secondary electron emission 
from the probe is neglected and the electron saturation current I^s is given by I^s = 
Aeen{l/A) -y/ (S/cb^g) / (Trme), using the random thermal current density. A^ and Ai specify 
the probe collecting areas for electrons and ions, respectively. Depending on the magnetic 
field strength, these areas can be differing for the two species, as stated in At any rate, 
to determine density and plasma potential from the measured quantities, electron and ion 
temperatures have to be taken into account, although Ti is frequently assumed to be equal 
to Tc. This is mainly due to the fact that in the edge region of fusion devices there is often 
no data available for the ion temperature. 

The measurement of electron temperature fluctuations by means of classical Langmuir 
probes requires a sweeping of a preferably complete probe characteristic. This results in a 
lower time resolution of the respective time series compared to data acquired by floating 
or negatively biased lis probe pins, although the method has undergone further develop- 
ment towards fast sweeping probes (jsl-jsl). Alternatively, triple probes (joj, 10|) or the 



harmonics technique ( lll-llSi) can be used to measure fluctuations of Tp. In addition, more 

n n 

sophisticated probes such as emissive probes in contrast to conventional cold 

probes, or ball-pen probes (jl5|, 16|) have been developed. These kind of probes are aimed 



at measuring the plasma potential directly and a combination of cold probes and ball-pen 
or emissive probes also allows a derivation of the electron temperature [l^]. Nevertheless 
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most probe measurements in the edge of large fusion devices are still based upon data mea- 
sured by classical Langmuir probes. For estimations of the radial particle flux T,- = nvj-, the 
radial velocity is commonly calculated from gradients of the floating potential instead of the 
plasma potential 18| and the density is calculated from I^g using only average values for the 
temperatures. 

The aim of this paper is to use numerical gyrofiuid simulations to compare time series of 
density n, plasma potential $ and temperatures Te, T; with simulated values of ion saturation 
current lis and floating potential Vfl. The difference between the two corresponding datasets 
of $ and Vfj is to some extent comparable to the difference between emissive and conventional 
probe measurements, although the potential measured by the former still deviates from the 
actual plasma potential by a certain, albeit smaller temperature-dependent factor 1^ . 

For the simulations the nonlinear three-dimensional electromagnetic gyrofiuid turbulence 
code GEMR has been used, which comprises a six-moment gyrofiuid model for electrons 
and ions in a circular toroidal geometry and features energetic consistency ( 19|-j23|). The 
coordinate system in use consists of a flux surface label (x) defining the radial position, a 
field line label within the flux surface (y) and a position along the field line (s) 2^. For 
diagnostic purposes, the code delivers time series of fluctuating electron and ion densities, 
plasma potential, temperatures and parallel velocities, amongst others. Thus, the knowledge 
of these quantities including electron and ion temperatures allows the calculation of ion 
saturation current and floating potential. So the significance of temperature fluctuations 
with regard to density and potential measurements can be investigated. Moreover, the 
analysis can be performed at different simulated probe positions in the radial computation 
domain, which in our nominal case is (r/a = 1 ± 0.06), with a being the minor radius. 

Typical ASDEX Upgrade edge values have been chosen as input parameters for the 
simulation. The model can be regarded as global in the sense that there is a global variation 
of profiles, although the parameters are constant [l^. Although no turbulence code can 
as yet self-consistently achieve an H-mode, magneto-hydrodynamic ideal ballooning modes 
(IBMs), which are commonly assumed to cause edge localised modes of type I, can be 
simulated by incorporating ex per imental H-mode density and temperature pedestal profiles 
n(r) and T{r) as initial state 2l|. Ideal ballooning (ELM-like) blowouts in experiment and 



simulation are always connected with large fluctuations in density, temperature and plasma 
potential, so a comprehensive analysis of the temperature influence on floating potential and 
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ion saturation current in this case is a reasonable addition to investigations of simulations 
in saturated L-mode state. 

In the following sections the details of the comparison between plasma density and ion 
saturation current, plasma potential and floating potential and characteristics of the par- 
ticle flux in L-mode situation (section 2) as well as IBM blowout situation (section 3) are 
presented. Possible reasons for the varying discrepancy between the actual quantities and 
the simulated Langmuir probe measurements are also addressed and the role of temperature 
fluctuations is discussed. 

II. SATURATED L-MODE SITUATION 

The first situation corresponds to an operation in saturated L-mode. That means, the 
electron dynamical plasma beta /3e = (/WoPc) / is chosen low enough not to be ideal balloon- 
ing unstable (/3c ~ 9.4- 10~^) and electron and ion heat sources as well as density sources are 
set to a moderate level. The background mid-pedestal parameters for this case are Tg = 150 
eV and Ti = 180 eV for the temperatures, ne = n\ = 1.25-10^^ m~^ for the densities, B = 2.0 
T for the background magnetic field, Lt = L± = 3.0 cm for the perpendicular temperature 
gradient length and = 6.0 cm for the density gradient length. For the ion mass, the 
deuterium mass mo = 3670mc is used. Major torus radius and aspect ratio (major divided 
by minor radius) conform to ASDEX Upgrade values with R = 1.65 m and R/a = 3.3, while 
a circular flux-surface geometry is used. 

The radial domain, whose direction is defined as the x-direction in the code, is divided 
into Tlx = 64 grid points. These are not equally spaced in terms of radius but in terms 
of volume, which for the present purpose makes a slight but not decisive difference. The 
y-direction is divided into Uy = 512 grid points, the s-direction into = 16 grid points and 
the averaged grid size perpendicular to the magnetic field is 1.06ps x 1.39ps (with ~ 0.88 
mm). The analysed L-mode time series have a length of 11000 data points and a temporal 
resolution of 0.0826 /is (total duration 0.908ms). 

Especially in the outer scrape-off layer region of the radial domain (x ~ 33 — 64) the 
numerical (Arakawa)-scheme occasionally delivers unphysical negative density and temper- 
ature values in the presence of steep propagating gradients, arising from Gibbs oscillations. 
Therefore we added the absolute value of the largest negative value multiplied by an off- 
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FIG. 1: (a) Time series of the electron temperature T^^^. (b) Mean values and standard deviations 
of electron and ion density and n\ (black), plasma potential $ (blue), electron temperature 
Tgii (red) and ion temperature Tiy (green), normalised to the respective mean at x = 32, for the 
L-mode case at different radial positions. The simulation grid point x = 16 (black bold line in (a)) 
corresponds to the radial distance of ~ —14.8 mm w — 16.7/9s, x = 32 (blue line in (a)) to ~ 0.4 
mm ~ 0.5ps and x = 48 (red thin line in (a)) to ~ 15.2 mm ~ 17.2ps, measured from the separatrix 

set factor Aqs to all values of the time series at all positions, e.g. for the density, n^^x = 
ni^x + Aos|Min(ni)|, where fii^x is the uncorrected time series at x and the minimum of all 
values in the radial domain is used. As offset factor, Aos = 1-01 has been chosen. 

The simulation run exhibits an ion temperature gradient (ITG) driven interchange over- 
shoot quite at the beginning and then saturates rather quickly to an L-mode-like state. For 
the analysis and the plots presented in this section only the well saturated part after the 
initial transient has been taken into account. Fig. [T] exemplary shows the time series of Te\\ 
as well as mean values and standard deviations of Ue and ni, $, Tc|| and Tiy at different radial 
positions (inside, near and outside the separatrix). 

Both here and in the following parts of the L-mode section, plots showing the temporal 
evolution of a quantity are based on data taken near the outboard midplane, whereas plots 
of time averaged data are composed from the individual results of different toroidal positions 
to provide better statistics. The radial position outside the separatrix (~ 1.5 cm, red thin 
lines) corresponds approximately to the region of Langmuir probe measurements in ASDEX 
Upgrade, which are typically positioned few centimeters outside the last closed flux surface. 
A separate treatment of parallel and perpendicular temperatures (Te||, Tiy and Te±, Ti±) 
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arises from the construction of moments in the gyrofluid equations 25|]. However, in the 
present case there is no major difference between both components and only time series of 
the parallel temperatures have been used to calculate the quantities of the synthetic probe. 



A. Plasma density vs. ion saturation current 

The ion saturation current lis has been calculated from ion density and the temperatures 
of electrons and ions according to eq. [T] (with A^, = defining the cross section of the 
probe). As the characteristics of the n[, T^w and T[\\ time series are rather similar, there 
are no striking differences between the ion saturation current, which is essentially a product 
between them, with the sum of temperatures appearing as a square root, and the underlying 
quantities and it still shows a qualitatively comparable temporal evolution. 

For evaluations of the radial particle flux from experimentally measured probe data it 
is necessary to know, amongst others, the particle density. Although the conversion can 
easily be done using eq. [H in many experimental cases only average temperature values 
are available (e. g. evaluated by sweeping the Langmuir probe). Hence it makes sense to 
recalculate Ui from /is using averaged values of the simulated Tc|| and T[\\ time series for 
each x-position in the radial domain and to compare the results with the precise ri; time 
series. Both signals are plotted in fig. |2]for a position within the SOL. Inside and near the 
separatrix similar results have been obtained. The difference between nf"^^ and ni, given by 

is particularly large at points of maxima and minima, but rather small elsewhere. As n'^^^ = 
rii^ (Te + 71)/ (Tl"^^ + T^^^\ the ratio of the actual temperature values to their average values 
is the decisive factor. The average temperatures are smaller than the maxima, which results 
in larger values of n^^^ at these positions (the temperatures appear in the denominator in 
eq. [3]). For the minima, the same statement holds vice versa. 

B. Plasma potential vs. floating potential 

Fig. [3] shows a comparison of the plasma potential $ and the floating potential V-a for 
a position within SOL, computed by means of eq. [2l A clear difference is evident, as is to 
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FIG. 2: (a) L-Mode time series of code output n\ (red bold line) and n^'^^ (black thin line) calculated 
from lis using averaged values for T^n and Tiy . The data is plotted for a radial position ~ 15.2 mm 
w 17.2ps outside the separatrix. (b) Detailed evolution of a 160 /US time frame 

be expected, both in mean values and fluctuations. Comparing different radial positions, 
the mean offset between $ and Va is decreasing from the inside to the outside. In terms of 
mathematics, the difference is caused by a product of the electron temperature in energy 
units, divided by e, and a dimensionless quantity, which depends on the temperature as well: 



A 



real 



T 



A. 



(4) 



e / \ Ai\ Tc + T^y 2nmc I \ e 
If the temperature is assumed to be constant, i.e. all temperature fluctuations are ne- 
glected, the latter can be understood as average difference A^^^ between $ and Va (also 
referred to as a), normalised by e/{k-BTc) (see It can roughly be estimated for our nom- 
inal case with Tj = 1.2Te and = to be A^^^ ^ 2.79. Thus on average the normalised 
difference for a certain relation of Tc and Ti is constant and deviations are mainly due to 
temperature fluctuations. On the other hand, the background mid-pedestal temperature 
values used here are only reference values. For all radial positions, but especially within the 
scrape-off layer the simulated temperature relation can differ considerably from this. Con- 



sequently, A^^^ is subject to radial variations (see fig. 5(a)). Taking into account that the 
temperature is not constant results in a fluctuating time series A, which is plotted for one 



exemplary position within the SOL in fig. 4(a) The radial profile of its standard deviation 



is shown as light grey area in fig. 5(a) In the same figure, source and sink regions of the 
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FIG. 3: (a) L-Mode time series of $ (red bold line) and Vn (black thin line) for a radial position 
~ 15.2 mm ~ 17.2/3s outside the separatrix. (b) Detailed evolution of a 160 fis time frame 
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FIG. 4: (a) Normalised difference A between Boating and plasma potential for a position ~ 15.2 
mm ~ 17.2ps outside the separatrix, calculated from the Euctuating temperature time series (black 
solid line) and from constant temperature values averaged at the corresponding radial position 
(red dashed line), (b) Detailed evolution 130 /US of the actual difference A'''^^' between Boating and 
plasma potential (black line) and {(kBTe) / e) A^^^ (blue line) within the SOL 

GEMR simulation model are indicated by dark grey shadings. The radial positions affected 
by these mechanisms are removed in all subsequent profile plots, as they should be excluded 
from the interpretation. 

It is important to note, that A is the normalised difference between $ and Vq. It can 
be a useful quantity for rough evaluations, if only average temperatures are available. A 
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FIG. 5: (a) Radial profile of A'^^s (^l)lack line with circles), standard deviation of exact A time series 
(light grey area) and source/sink regions of the GEMR simulation model (dark grey shading), (h) 
Radial profile of (A'''^^^) (black line with circles), (Tc) (red line with squares) and the respective 
standard deviations (light grey area). The location of the separatrix is indicated by the dotted line 

quantitative estimation of the actual difference requires the calculation of A multiplied by 
{kBTc)/e (eq. H]) for every time step and for all positions. Thus the shape and phase of the 
temperature time series, which is varying depending on the measurement position, has a 
direct influence on the difference between plasma and floating potential. 

As the mean temperature (Tc) strongly decreases within the radial simulation domain, 
the actual difference A^^^^ (oc TcA) between floating and plasma potential is larger inside the 



separatrix and becomes smaller in the SOL (fig. 5(b) ). The difference in terms of fluctuation 
amplitudes seems to be somewhat larger inside and near the separatrix, which is due to large 
temperature fluctuations in this region, but a significant distinction of the general shape is 
evident at all radial positions. 

In order to emphasise the role of the electron temperature fluctuations, a comparison 
between the temporal evolution of A'^'^^' = ((/cBTe)/e)A and the quantity {{kBTe)/e)A^^^ is 



shown in fig. 4(b) There are only faint differences visible, indicating the predominance of the 
fluctuations of Tc compared to the fluctuations of A and, consequently, to the fluctuations 
of Ti. This behaviour can be observed in the entire radial domain, with a slightly decreasing 



difference between the two terms plotted in fig. 4(b) from the inside to the outside and 
therefore a decreasin g im portance of a precise A. 



As stated in ref. 



171], the differences in fluctuations of Vg and $ are determined not 
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FIG. 6: (a) Radial proGle of [RMS(Vfi)]2 - [RMS($)]2 (black line with circles), which is equal to 
the difference of [RMS(A''°^')]2 (blue solid line) and 2($A'^^i) (red line with squares) The tilde 
indicates fluctuating parts and RMS root mean square values 

only by fluctuations of the electron temperature, but also by the phase relation between 
temperature and potential. This becomes evident from the difference of the corresponding 
root mean square values, 

[RMS(Vfl)]2 - [RMS(i>)]2 = [RMS(A''*^'^')]^ - 2(4>A''^^^). (5) 

The tilde indicates fluctuating parts and the RMS of, for example, $ is given by ($^)^/^, 
which in the case of fluctuations with zero mean reflects the standard deviation (except for a 
slightly different normalisation). Fig. IH] shows radial profiles of the left-hand side of eq. Oand 
the two terms on the right-hand side. It is apparent, that although the contribution of A'''^^^ 
(oc Te) is clearly dominant, the influence of the cross phase between $ and A'''^'^^ represented 
by ($A''°^^), is of non- negligible magnitude and radially varying, with largest absolute values 
near the separatrix and within the SOL, but rather small values inside the separatrix. A 
comprehensive analysis of the phase relation between temperature and potential at different 
radial positions by means of time-resolved Wavelet methods is currently in progress and will 
be presented in a future work. 

The present calculation of the floating potential using simulated temperature time series 
depends to a certain degree on the artificial offset factor Aqs, which has a decisive influence 
on the temperature averages and therefore also on the mean offset between $ and Va. Apart 
from that, time series of A are affected, as the logarithm in the conversion equation eq. [2] 
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reacts quite sensitive on variations of the temperature minima, whose difference to zero is 



defined by Aos- As a consequence, not only the distinct negative peaks in fig. 4(a), but 
also the radial profiles of A'^^s and Stddev(A) change their characteristics depending on 
the chosen offset factor. Since Aqs has no effect on temperature fluctuations, its impact on 
the fluctuating part and the standard deviation of A"^"^' and on the comparison in fig. 4(b) 
is weak. That implies, that although the shape of the A time series and the mean offset 
between $ and is affected by this numerical issue and the artificial offset, the influence 
of electron and ion temperature fluctuations on calculations of the floating potential can 
nevertheless be studied. This applies even more to calculations of the radial particle flux, 
because only potential gradients are of importance there. A prerequisite for any consider- 
ations about the influence of small fluctuations on potential measurements is the question, 
whether an experimental probe is in principle able to react almost instantaneously to tem- 
perature fluctuations. This can indeed be assumed, because density fluctuations, which have 
similar characteristics and a comparable ratio of fluctuations to mean, can be detected. 

C. Radial particle flux 

One of the most important quantities for the statistical analysis of experimental mea- 



surements is the turbulent fluctuation-induced averaged radial particle flux ( 26|- 29| and 
Q), defined by 

Tr = (^^^^r) ~ — ^ , (6) 

where the fluctuating radial drift velocity is assumed to be the radial component of the 
fluctuating E x B velocity, Vr ~ "^^exb ~ Epoi/B. Appropriate averaging is indicated by 
(■). The poloidal electric field can be approximated by taking the difference between two 
potential measurements divided by their poloidal separation distance. In our simulations, 
this distance is given by di2 ~ 2.5 mm. Replacing the density by the expression for the ion 
saturation current lis from eq. [1] yields 

^'"^si {jTe]! k^{Tl + Tif') ( '''du ) ■ 
Both, the synthetically measured particle flux (using Va and densities inferred from 
lis and averaged temperatures) and the real particle flux FJ!®^' (using code outputs $ and rig) 
have been calculated and their temporal evolution, representing the instantaneous particle 
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FIG. 7: Radial particle Eux Tj- calculated from plasma potential and density (Tl'^'^^, red bold line) 
and from Boating potential and ion saturation current (T™, black thin line). The data is plotted 
for a radial position ~ 15.2 mm w 17.2ps outside the separatrix. (b) Detailed evolution of a 160 //s 
time frame. 
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FIG. 8: (a) Radial profiles of the time averaged particle flux (Fr) (averaged across the entire 
simulation time), (b) standard deviation Stddev(ri.). Real values are plotted as red line with 
squares, measured values as black line with circles 

flux, is plotted in fig. [71 The time series of F™ exhibit much larger fluctuations than that 
of FJ!"^' at all radial positions, but especially inside and near the separatrix, which is clearly 



visible in the radial profile of the standard deviation of the time-averaged flux (fig. 8(b) ). On 
average, Stddev(F™) is ~ 7 times larger than Stddev(FJ!'^'^^). Magnitudes and radial profiles 
of the mean fluxes (FJ^"^') and (F™) depend on the length of the investigated time range, 
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but significant differences can be observed in any case (fig. 8(a), (F™) ~ 3.5 times larger 
than (rj!'^'^')). The use of an averaged temperature value in the density calculation only gives 
rise to small deviations (see fig. E]) and is not responsible for the large discrepancy in terms 
of fluctuations. Hence it must primarily be due to the radial velocity calculated from the 
gradient of the floating potential. 

The equivalent use of the floating potential instead of the plasma potential, which is the 
usual procedure to estimate the particle flux from Langmuir probe measurements, is based 
on the assumption, that the difference between plasma and floating potential is roughly 
constant at two closely spaced locations. As only the difference between two adjacent time 
series is relevant, the error caused by using instead of $ is expected to be small, depending 
on different deviations from this constant at the two locations. Such deviations have been 
discussed in the previous section, but for radially varying positions, whereas here different 
positions in the y-direction of the simulation model are of importance. Fluctuations of time 
series at these positions correspond to fluctuations of two poloidally separated positions, 
because ^ k_i_. For Va = $ — [{kBTc)/e]A, the difference between approximations of 
plasma potential gradient and floating potential gradient is given by 



$1 - $2 14,1 - ■^4,2 1 



di2 di2 di2 

= (Te,iAi-Te,2A2). (8) 

"12 e 

This expression vanishes, if Te^i = Te,2 and Ai = A2 is assumed for the two adjacent 
positions. According to fig. IH which shows 5$ = $1 — $2 and 6Vq = Va^i — V{i^2, this is 
clearly not the case, so there must be a significant difference between these quantities. The 
second part of equation eq. |8] can be further expanded by splitting the respective quantities 
into mean values and fluctuations: 

(Te,iAi - Te,2A2) = [((Te,i)(Ai) - (Te,2)(A2)) + ((Te,i)Ai - (Te,2)A2) 

+ ((Ai)Te,l - (A2)Te,2) + (Te.iAi - Te,2A2)] . (9) 

A small offset between the two time series of 5$ and 5Vfl is caused by slightly different 



mean values of Tg^i and Te,2 as well as Ai and A2 (fig. 10(a) ), which is expressed by the first 
term on the right-hand side of eq. [91 The remaining terms, the second last of which providing 
by far the largest contribution, are related to the strong fluctuations of 6Vf[. These originate 
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FIG. 9: (a) Difference between two time series $i and $2 fred bold line) and between Vq^i and 
Vfi^2 (black thin line), collected at the poloidally adjacent positions used for the estimation of Vr- 
The data is plotted for a radial position ~ 15.2 mm w 17.2/9g outside the separatrix. (b) Detailed 
evolution of a 160 /US time frame. 



from differences in the fluctuating parts of Te^i and Te,2 (fig- 10(b) ) and, accordingly but of 
less importance, Ai and A2. 

To point this out, and the differences 5$ and 6Vf{ of spatial adjacent time series have 
been calculated for the artificial case of temperature time series, which have equal mean 
and therefore no offset at the two positions, and whose difference in fluctuations is the same 



as before (fig. 10(b)), but with much smaller amplitude (1/10). The result is plotted in 



fig. 10(c), where 5V{i is in substantial better agreement with 5$ than in fig. 9(a 



In order to investigate the impact of temperature fluctuations on statistical properties, 
probability distribution functions (PDF) have been computed. As shown in fig. [TT] and 
fig- [121 there is a different behaviour inside and outside the SOL. Inside and near the sepa- 
ratrix, the PDF of the real particle flux features a larger skewness at most radial positions, 
whereas the kurtosis is largely of about the same magnitude. In the scrape-off layer, from a 
radial position of ~ 6 mm ^ Q.Sps outside the separatrix, both the skewness and the kurtosis 
of ri'^'^^ are smaller than that of F™, which indicates, that in the time series of F™ strong 
deviations from the mean value occur even more frequently than in FJ®*^^. 

Above investigations show, that both flux time series and their statistics are differing to 
a significant degree depending on whether temperature fluctuations are taken into account 
for the calculations. This seems to be the case not only for measurement positions inside 
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FIG. 10: (a) Temperature time series T^^i (black line) and T^^i (blue line), collected at the two 
poloidally adjacent positions, (b) Difference between the Euctuating parts of the two temperature 
time series, (c) Difference between the two spatially adjacent time series for ^ (red bold line) and 
Vfi (black thin line) for the artificial case of small differences of the Euctuating temperature time 
series (see text). All plots show the detailed evolution of a 160 /is time frame for a radial position 
~ 15.2 mm ~ n.2ps outside the separatrix 

and near the separatrix, but also within the SOL, where Langmuir probe measurements are 
typically located. Although the quantitative results of our numerical simulation might differ 
from the experimental situation in a fusion device and should not be taken as a reference, 
the importance of the fluctuating electron temperature is nevertheless evident. It should be 
mentioned, that variations of the offset factor Aqs, which is used to compensate numerical 
errors of the simulated time series, only have small influence on this and do not change the 
overall flux features. 



III. IDEAL BALLOONING MODE BLOWOUT 

As a next step, a simulated ELM type-I like ideal ballooning mode (IBM) situation 
has been analysed, which exhibits a large interchange blowout. A sudden burst connected 
with enhanced fluctuations is clearly visible in the time series of Teii, which is shown as 



an example in fig. 13(a), but of course also in the signals of rie and rii, $, and Tiy. In 
the aftermath of the blowout the standard deviations of the temperatures and densities are 
lower compared to the L-mode case at most radial positions, but the standard deviations of 
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FIG. 11: PDF of the radial particle Eux Fr, calculated from plasma potential and density (red bold 
line) and from Boating potential and ion saturation current (black thin line). The data is plotted 
at the radial positions (a) inside, (b) near and (c) outside the separatrix as specified in iig. [I] 
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FIG. 12: Radial profiles of (a) the Skewness and (b) the Kurtosis of real (red line with squares) 
and measured particle Eux F^ (black line with circles) 



the potential time series are larger (fig. 13(b)). The background temperature and density 
parameters are Te = 300 eV, Ti = 360 eV and ric = n\ = 2.5 ■ 10^^ m~'^. Gradient lengths 
and background magnetic field remain unchanged compared to the L-mode case of the 
previous section. This yields /3e ~ 4.0 ■ 10~^ and allows the occurrence of an IBM instability, 
according to the MHD ideal ballooning criterion am = > s (with /3 = 2/3c, safety 

factor q = 1.5 + 3.5(r/a)^ and magnetic shear parameter s = {r / q){dq/ dr) — ?■ Sa = 1.4 at 
the separatrix). Heat and density sources are set to a small value, in order to provide an 
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FIG. 13: (a) Time series of the electron temperature T^n . (b) Mean values and standard deviations 
of electron and ion density and n\ (black), plasma potential $ (blue), electron temperature 
TqIi (red) and ion temperature Tj|| (green), normalised to the respective mean at x = 32, in the 
aftermath of an IBM blowout at different radial positions. The simulation grid point x = 16 (black 
bold line in (a)) corresponds to the radial distance of ~ —14.8 mm ~ —ll.Sps, x = 32 (blue line 
in (a)) to ~ 0.4 mm ~ 0.34/3s and x = 48 (red thin line in (a)) to ~ 15.2 mm ~ 12.2ps, measured 
from the separatrix 

IBM blowout as pure as possible and to prevent distortions caused by other effects. The 
averaged grid size perpendicular to the magnetic field is 0.75ps x 0.98ps (with ^ 1.25 
mm). The analysed IBM time series have a length of 15000 data points and a temporal 
resolution of ~ 0.0167/is (total duration 0.25ms). All signals of this section are acquired 
near the outboard midplane. Further details on IBM simulations with GEMR can be found 
in ref. 21 1. 



A. Plasma density vs. ion saturation current and plasma potential vs. floating 
potential 

Just as in the L-Mode case discussed above, there are no major differences between the 
time series of density and ion saturation current and both feature similar characteristics. If 



the density is inferred from Jjs by means of averaged temperature values (fig. 14(a)), the 
main features are preserved, apart from the time- dependent offset mentioned in subsection 
III At which results from averaging. As the use of time-averaged values is more defective for 
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large density and temperature variations, it is most pronounced during the blowout. 

$, Vfi and A^^^^, the actual difference between plasma and ffoating potential, have also 



been calculated for the IBM case (fig. 14(b) and fig. 15(a)). An individual investigation of 
A has been omitted, since according to section IllBI the inclusion of the factor plays an 
important role. The predominance of Tg compared to A is still present, as a comparison 



between A''°'^' and ((/cB^c)/e)A'^^^, shown in fig. 4(b) for the L-mode situation, yields sim- 
ilar results for IBM time series. Differences between fluctuations of $ and Va are clearly 
noticeable at all radial positions, but due to the limited lengths of the IBM time series no 
reliable results for the standard deviation could be obtained, in particular before and during 
the blowout. After the blowout, the profile of the standard deviation seems to be similar to 
the L-mode case. Inside the separatrix, the IBM blowout leads to a strong reduction of the 
temperature, whereas near and beyond the separatrix it is increased. This is reflected in the 



temporal evolution of A'''^'^^. Fig. 15(b) shows radial profiles of (A'''^'^^) and (Tc), separately 
averaged across the time range immediately before and after the blowout. Inside the sepa- 
ratrix the time average of the actual difference and the temperature is considerably larger 
before the blowout than afterwards. From a point near but not across the separatrix it is just 
the other way around. The decrease or increase of A'''^^' after the blowout has implications 
if the plasma potential is calculated from Va using temperatures averaged across the entire 
time domain, as the resulting time series will be either underestimated before the blowout 
and overestimated afterwards or vice-versa. 



B. Radial particle flux 

Also in the IBM case, the instantaneous "measured" radial particle flux F™ calculated 
from /is and Va shows much larger fluctuations than the real flux Tl'^^^. This is again primarily 
due to the use of Va instead of $. The averaged temperature in the calculation of density 
from /is leads to a further amplification of the fluctuation amplitudes. The magnitude of 
strong flux peaks during the blowout, which varies considerably at different radial positions, 
has a decisive impact on radial profiles of mean fluxes and standard deviations, if these are 
averaged across the entire simulation time. Hence the time spans during and after the IBM 
blowout should be treated individually. In fig. [161 the time range 29.1 fis — 97.3 /is has been 
chosen, which includes the large peaks of the blowout at all radial positions. Fig. [16] shows 
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FIG. 14: (a) IBM time series of code output n[ (red bold line) and n^^^ (black thin line) calculated 
from lis using averaged values for and Tiy . (b) IBM time series of $ (red bold line) and Vq (black 
thin line). The data is plotted for a radial position ~ 15.2 mm ~ 12.2/3s outside the separatrix. 
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FIG. 15: (a) Actual difference A^'^'^^ between Boating and plasma potential for the radial positions 
~ 14.8 mm ~ ll.Sps inside (black line) and ~ 15.2 mm ~ 12.2ps outside the separatrix (blue line), 
(b) Radial profile of (A''®^') before (black line with circles) and after the blowout (red line with 
squares). Profiles of (Tg) are plotted as dotted line. 

profiles from the time range 97.3 /is — 239.5 /is. Small displacements of the transition point 
yield slightly different profiles, but the general characteristics remain unchanged. 

The larger amplitudes of F™ compared to Tl'^^^ can be observed at most positions in the 
radial profiles of mean values and standard deviations. The ratios between measured and 
real values are approximately comparable to the values given in the L-mode section. As 
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FIG. 16: (a) Radial profiles of the time averaged particle flux (T^), (h) standard deviation 

Stddev(rr) during tlie IBM blowout (29.1^8 97.3^sj. Real values are plotted as red line 

with squares, measured values as black line with circles 




30 40 
Radial position x 

(a) 



30 40 
Radial position x 

(b) 



FIG. 17: (a) Radial profiles of the time averaged particle Eux (Fr), (b) standard deviation 

Stddev(Fr) after the IBM blowout ('QT.S^s 239.5/isj. Real values are plotted as red line with 

squares, measured values as black line with circles 



shown in fig. p!6(a)j {Tf") and (F;;^^!) differ at some radial positions not only in absolute 



values but also in sign, which is a significant distortion of the actual situation. However, 
the respective mean values are based only on a relatively small time interval and should not 
be overemphasised. Interestingly, the fiuctuation amplitudes and the standard deviation of 
F™ and Tl"^^^ approach a local minimum near the separatrix. This is especially pronounced 



during the blowout (fig. 16(b)), but is still present afterwards (fig. 17(b)) 
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As our simulation code can only provide data for the temporal evolution of one single 
IBM blowout and its aftermath, the fluctuations covered by the time series are not sufficient 
to allow a more detailed statistical analysis including probability density functions and 
statistical moments. 

IV. CONCLUSIONS 

From the above investigations can be concluded, that neglecting any temperature fluctu- 
ations in our computations of probe measurements results in signiflcant differences between 
the synthetically measured values and the actual quantities at all radial positions of the 
simulation domain, both in a saturated L-mode situation and in simulations of an IBM 
blowout. This holds especially for the floating potential and its use in calculations of the 
fluctuation induced particle flux. The latter shows a considerable discrepancy not only in the 
temporal evolution but also in statistical properties, if spatial variations of the temperature 
fluctuations are not taken into account. However, fluctuations of the electron temperature 
seem to be of much greater importance than ion temperature fluctuations. 

Compared to the saturated L-mode situation, the investigation of time series involving 
an IBM blowout did not yield a major alteration of the ratio between real and measured 
quantities, apart from unsurprising changes due to the large peaks. 

Although a realistic implementation of a virtual Langmuir probe in terms of geometry 
and particle processes would require the use of a kinetic model, the results from our gyrofluid 
simulations can nevertheless be regarded as relevant for experimental measurements, as all 
comparisons are performed within the simulation. In this respect they may serve as an 
indication to the particular importance of electron temperature fluctuations, whose neglect 
in evaluations of experimental data might lead to substantial deviations from the actual 
quantities. 
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